Short time dynamics of a two-dimensional majority vote model 
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Short time Monte Carlo methods are used to study the nonequilibrium ferromagnetic phase 
transition in a majority vote model in two dimensions. The existance of an initial critical slip 
regime is verified. The measured values of dynamic exponents z = 2.170(5) and 6 — 0.191(2) are in 
excellent agreement with those of the kinetic Ising model universality class. 
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The identification and caracterization of universality 
classes for nonequilibrium systems are far less settled 
than in the case of systems in thermal equilibrium. Nev- 
ertheless, many model systems with microscopic irre- 
versible dynamic rules (no detailed balance) and two 
states per site have been found to fall in the Ising model 
universality class, as far as the static behaviour is con- 
cerned [1-8]. Grinstein, Jayaprakash and Yu He 1^ have 
argued that, provided the rules are up-down symmetric, 
both the statics and the (long time) dynamics of those 
models are the same as those of the kinetic Ising model. 
We are however not aware of any direct determination of 
the dynamic z-exponent for such models. In the present 
work we have applied the recently proposed early time 
dynamic Monte Carlo technique jlOj, |11|, to inves- 
tigate the dynamic behaviour of a model in the above 
mentioned class - the majority vote model (whose rates 
can be seen as a combination of two Glauber dynam- 
ics in contact with two heat baths at different tempera- 
tures) ijl^ , [Q . Previous studies of short time dynamics 
were concerned either with equilibrium systems, like Ising 
0, H or Potts |17) models, or with a 

nonequilibrium phase transition in a distinct universality 
class 

Janssen, Schaub and Schmittmann have shown 
that, when a system with relaxational dynamics is 
quenched from T ^ Tc to Tc, the early times of evo- 
lution also display universal behaviour. A new indepen- 
dent exponent, 9, associated with the anomalous dimen- 
sion of the initial order parameter, was introduced to 
describe the system behaviour during this critical initial 
slip regime. Denoting by mp the initial magnetization 
(0 < mo <C 1), this regime is found in the time range 
where tmic is some microscopic time 
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< t < m. 
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and xo — 9z + (3 /v {(3 and v are the equilibrium critical 
indices). The magnetization (m(t) = N^^ < di >) 
increases with time as 



m(<) moi^ 



(1) 



is also related to the decay of the auto-correlation func- 
tion from a disordered initial state, 

A{t)^t-^ (2) 

with A = d/ z — 9 in d space dimensions. The relation 



between short time dynamics and damage spreading was 
recently clarified by Grassberger p9|]. 

For a finite system it is expected |10| [pd| that the mo- 
ments of the order parameter, m*-'^-' {kth moment of the 
magnetization), have the scaling form, 

m^^) (i, T, L, mo) - h-^^'^m^''^ (fo^^i, b-^'^T, bL, fo'^^mo) 

(3) 

where t — (T ~ Tc)/Tc and the initial correlation length 
is null. Setting r = and the arbitrary scaling factor 
b ~ i^/^, one obtains from eq.(3). 



m(t,L,mo) - t'^/^'Fit/tL^t/to) 



(4) 



where t^ ^ and to = ' °. Following the scaling 
relations for the magnetization and its higher moments, 
it is possible to infer that the time dependent Binder 
cumulant |20|, 



U{t,L) = 1 



7,(4) 
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obeys. 



Uit,Li) = U{b~H,L2) 



(5) 
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for T — tuq = and two system sizes (Li and L2) with 
b = L2/L1. The exponent z can be obtained from a data 
collapse with a time rescaling factor b~^. Since only early 
times are considered, this is a rather efficient method to 
extract z. Once z is known, the static exponent P/i^ is 
recovered from a similar scaling analysis of m*^^) 

Starting with random initial configurations (with mo = 
0) and following the evolution at Tc of the spin auto- 
correlation function, 
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the power-law decay, eq.(2), is visible after a short tran- 
sient regime and A (and therefore, 9) can be obtained 
p6[ . A direct measurement of 9 is possible from eq.(l) 
- the samples are then prepared with a sharply defined 
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small value of mo. After a few Monte Carlo steps (MCS) 
a straight line appears in the log-log plots and 9 is com- 
puted from its slope. 

The 2-state isotropic majority vote model is defined 
by a set of 'voters' or 'spin' variables {ci} taking the 
values +1 or —1 and evolving in time by a single spin- flip 
like dynamics with a probability Wi given by 



l-a^il-2q)S iJ2<^^^ 



(8) 



where S{x) = sign{x) if x ^ 0, S{x) = if a; = and 
the sum is over nearest neighbours of ct^. The control 
parameter q plays the role of temperature in equilib- 
rium systems and measures the probability of aligning 
antiparallel to the majority of neighbours. In two dimen- 
sions this model has a ferromagnetic stationary phase for 
< q < (Zc undergoing a second order phase transition to 
a paramagnetic phase at qc {qc — 0.075(1) for a square 
lattice |i|L B). The static critical behaviour is Onsager- 
like 0, [g. According to Grinstein et al's arguments 

, its dynamic critical behaviour is the same as model A 
(Ising) |2l] and therefore Janssen et al's renormalization 
group analysis should also apply to such a nonequilib- 
rium model. The aim of this communication is to report 
a direct confirmation of the above conjecture from the 
results of a short time dynamics study of this model. 
The exponents z and 6 are found to be indistinguishable 
from the corresponding Ising values, z — 2.172(6) [ p^ , 
z = 2.1665(12) m and = 0.191(3) Hi and the exis- 
tance of an intermediate scaling regime is verified. 

Simulations were carried out for square lattices of side 
L = 16, 32, 64 and 128 with periodic boundary condi- 
tions. Random initial configurations with mo = were 
used in the study of U{t) and A{t), whereas a small ex- 
cess of plus spins, randomly distributed on the lattice, 
was taken to produce a selected value of toq 7^ 0. In 
order to prepare a sample with a precise magnetization 
and negligeable correlation length, we generated a lattice 
state with equal probability of occupation for both spin 
states and then flipped the spin at randomly chosen sites 
until the desired magnetization was obtained. The lat- 
tice was updated by flipping randomly picked spins ||2^ ] 
with probability given by eq.(8) with q = qc ~ 0.075. 
The evolution was followed for up to 1000 MCS. Aver- 
ages were performed over a large number of histories (up 
to 4 X 10^ independent initial configurations). 

In Figure 1 Binder's cumulant, U{t,L), is displayed 
against t/L^ for different values of system size and ini- 
tial magnetization mo = 0. The value of the z-exponent 
obtained from the best collapse was 2.170(5), which 
is in very good agreement with the ones obtained by 
Grassberger from a damage spreading study and by 
Nightingale et al from a new method (variance re- 
ducing) Monte Carlo algorithm for the two dimensional 
Ising model with nonconservative dynamics. Figure 2 
shows a plot of log m{t) against log t for mo = 0.03125 
and different system sizes {L — 16, 32, 64 and 128); for 



comparison a straight line with slope 0.191 is also drawn. 
The dependence of 9 on toq was analised (see Table I). 
A linear extrapolation to mo = yields = 0.191(2). It 
is clear from Figure 2 that 9 can be estimated from the 
study of very small system sizes. 
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FIG. 1. The figure shows a collapse plot of the Binder cu- 
mulant U{t,L) as a function of t/L" for L = 8(*), 16(D), 
32(o) and 64(x) with z = 2.170. The inset shows U{t,L) 
against t (same symbols). 



In our case, even for L = 16, the power-law behaviour 
lasts for two decades of MCS with an exponent close to 
our best value. This shows that for the measurement 
of 9 the finite size effects are not important, in con- 
trast with the behaviour of the auto-correlation function, 
where higher values of L are necessary. From Figure 2 we 
can also obtain an estimate of the crossover time tc when 
the magnetization changes to the decreasing power-law 
behaviour {f^^'^^), before entering the ultimate expo- 
nential regime. It is clear that for t < tc the magneti- 
zation presents a power-law increase for all values of L 
with no significative finite size effects. It is also remark- 
able that this power-law appears at a very small time 
(W ~ 1 MCS). 

In Figure 3 we show a double-log plot of the auto- 
correlation function, A(t), as a function of time for dif- 
ferent values of L. To get the critical exponent A we 
discarded the first 10 MCS. 

TABLE I. Exponent ^ as a function of mo for L — 32. The 
value obtained for mo is ^ = 0.191(2), in good agreement 
with the ones in literature 11 ql, llSl. 
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0.023437500 


0.1870(20) 


0.031250000 


0.1856(20) 


0.058593750 


0.1809(20) 


0.080078125 


0.1783(20) 


mo 


0.191(2) 
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FIG. 2. Log-log plot of the time evolution of the magne- 
tizatiom m{t) for system sizes L = 16, 32, 64 and 128 from 
bottom to top. The initial magnetization was mo — 0.03125. 
The dashed guide line has slope 0.191 and was plotted for 
comparison. 



-1.0 




_4.0 I ' ' ' ' ' 1 

1.0 2.0 3.0 4.0 

logt 

FIG. 3. Log-log plot of the auto-correlation function A{t) 
for system sizes L = 32, 64 and 128 (bottom to top) with ini- 
tial magnetization mo — 0. The dashed guide line has slope 
0.735. 



The best value obtained for A was 0.735(5), which 
agrees with recently published results [Q. The scaling 
relation A = d/z — is well obeyed for the values obtained 
before, within the statistical errors. 

In summary, we have investigated the dynamic be- 
haviour of a critical nonequilibrium model, the two di- 
mensional majority vote model, making use of the early 
time dynamic Monte Carlo method. By following the 
time evolution of the magnetization. Binder cumulant 
and time auto-correlation function for systems of vari- 
ous sizes, we were able to calculate numerically the val- 
ues of the exponents z = 2.170(5), 9 = 0.191(2) and 
A = 0.735(5). These values are in very good agreement 
with their corresponding two dimensional Ising results - 
a direct confirmation of the stability of the kinetic Ising 
fixed point with respect to irreversibility of the micro- 
scopic rates (of a certain kind). The effect of the ab- 
sence of detailed balance has probably to be sought in 
properties like the cluster structure and the dynamics of 
pattern formation p^ . Another dynamic exponent, 0i, 
the global persistent exponent, was recently introduced 
by Majumdar et al |2^. It measures the persistence of 
the sign of the magnetization and is related to 9 for a 
Markovian system | |25[ | . Evidence of non-markovianicity 
was reported for a nonequilibrium mode l ||2q| but the sit- 
uation is unclear in the Ising case ||27[. The persis- 
tence probability of the majority vote model is currently 
being investigated. 
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